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Abstract. We describe the results high resolution, hydrodynamic calculations of neutron star mergers. The 
model makes use of a new, nuclear equation of state, accounts for multi-flavour neutrino emission and solves the 
equations of hydrodynamics using the smoothed particle hydrodynamics method with more than 10 6 particles. 
The merger leaves behind a strongly differentially rotating central object of ~ 2.5 M together with a distribution 
of hot debris material. For the most realistic case of initial neutron star spins, no sign of a collapse to a black hole 
can be seen. We argue that the differential rotation stabilizes the central object for ~ 10 2 s and leads to superstrong 
magnetic fields. We find the neutrino emission from the hot debris around the freshly-formed, supermassive 
neutron star to be substantially lower than predicted previously. Therefore the annihilation of neutrino anti- 
neutrino pairs will have difficulties to power very energetic bursts (3> 10 49 erg). 



INTRODUCTION 

There is growing observational evidence that the sub- 
class of long Gamma Ray Bursts (GRBs) is related to 
star forming regions (e.g. [4]). While this connection is 
immediately evident for the short-lived progenitors of 
collapsars, the question of whether neutron star binaries 
merge close to star forming regions or not is not a set- 
tled one. Bloom et al. (2002), for example, argue against 
compact object mergers as central engines of (at least the 
subclass of the long) GRBs. Belczynski et al. (2001), 
however, claim to have identified new formation chan- 
nels that lead to classes of very tight, short-lived neutron 
star binaries. They find typical inspiral times of the or- 
der 10 6 years and therefore neutron star binaries would 
merge very close to their birth places, even if equipped 
with high systemic velocities due to kicks in asymmetric 
supernova explosions. 

Recent afterglow observations of long bursts suggest that 
they are beamed and require Ey ~ 5 -10 50 erg [8]. The en- 
ergy requirements for the subclass of short bursts are so 
far essentially unconstrained. Due to their gravitational 
binding energy of several times 10 53 erg, neutron star 
binaries certainly do possess the energy reservoirs nec- 
essary to power a (long) burst, however, how to trans- 
form the available energy into gamma rays is still far 
from being clear. Among the suggested mechanisms are 
magnetic energy extraction processes (e.g. [6, 11, 22]) 
and the annihilation of neutrino-antineutrino pairs emit- 
ted from the hot neutron star debris during the coales- 
cence [7]. 



HIGH RESOLUTION SIMULATIONS OF 
THE MERGER EVENT 

To study the possible role of neutron star coalescences 
for either long or short bursts we have performed de- 
tailed high-resolution simulations of the merger event 
[15]. To solve the equations of fluid dynamics we apply 
the smoothed particle hydrodynamics method (SPH) to- 
gether with a largely improved artificial viscosity tensor 
[14]. We treat the self-gravity of the neutron star fluid in 
a Newtonian way, but we add the forces emerging from 
the emission of gravitational waves to drive the system 
towards coalescence. The microphysical properties of the 
hot and dense neutron star matter are described using 
an equation of state (EOS) that is based on the tables 
of Shen et al. [19, 20]. We have extended the EOS to 
the low density regime via a gas consisting of neutrons, 
protons, alpha particles, electron-positron pairs and pho- 
tons. This new EOS covers the whole relevant parameter 
space in density, temperature and electron fraction. To 
allow for cooling and compositional changes, we have 
implemented a detailed neutrino treatment that accounts 
for three neutrino flavours (v^V,,, and v v , which is col- 
lectively used for the four heavy-lepton neutrinos) and 
takes the relevant emission processes (lepton capture on 
nucleons, electron-positron pair annihilation and plas- 
mon decay) into account. The opacities are calculated 
from the absorption processes of the electron-type neu- 
trinos on nucleons and scattering off nucleons and nuclei, 
the latter process becoming the dominant opacity source 
even for moderate mass fractions of heavy nuclei. The 
calculations are performed efficiently on shared-memory 
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FIGURE 1. Coalescence of a corotating binary system of 
1.4 Mq per star. Colour-coded is the matter density in the 
orbital plane. More than 10 6 SPH-particles were used for this 
calculation. 



parallel computers with more than 10 SPH-particles. 
For a more detailed description of the input physics and 
the computational methods we refer to Rosswog et al. 
[15, 16]. 



RESULTS 

In this latest set of calculations we explore systems 
with two initial spin configurations (for details see [15]): 
tidally locked, corotating systems for the ease of con- 
structing initial equilibrium models and systems without 
initial neutron star spins. The latter ones are the most 
relevant spin configurations since the intrinsic neutron 
star viscosity is too low to lead to a tidal locking during 
the short phase where the binary components undergo 
tidal interaction [3, 10]. We find central objects of 2.3 
to 2.6 M© which are strongly differentially rotating. For 
the most realistic, irrotational case the maximum den- 
sity does not even reach the initial density of a single, 
cold non-rotating neutron star. Apart from the thermal 
pressure and the differential rotation there may be fur- 
ther effects that, at least temporarily, stabilize the merger 
remnant against collapse to a black hole: the presence 
of non-leptonic, negative charges together with trapped 
neutrinos [13], for example, can substantially increase 
the maximum possible mass. One also expects magnetic 
seed fields to be amplified in the differentially rotat- 
ing remnant to enormous field strengths (~ 10 17 G) [6]. 
Fields of this strength can substantially modify the struc- 
ture of the central object and provide additional support 



against collapse [5]. Due to our ignorance of the high- 
density equation of state it cannot be ruled out that the 
end product of the coalescence is a stable supermassive 
neutron star of ~ 2.8 M . It seems, however, more likely 
that the central object is only temporarily stabilized and 
once the stabilizing effects weaken (e.g. neutrinos have 
diffused out after ~ 10 s, magnetic braking has damped 
out differential rotation) collapse to a black hole will set 
in. The time scale until collapse is difficult to determine, 
since it depends sensitively on poorly known physics and 
on the specific system parameters. We expect the most 
important effect to come from the differential rotation 
[12, 1] and therefore the collapse time scale to be set by 
the time it takes the remnant to reach uniform rotation. 
Assuming the dominant effect to come from magnetic 
dipole radiation (the viscous time scale is estimated to 
be ~ 10 9 s [18]), this time is given by 

18c 3 M 
Xc ~ 5fi2R 4 co 2 

where M,B and R are mass, magnetic field and radius 
of the central object. Such an object, an at least tem- 
porarily stabilized, supermassive neutron star with enor- 
mous magnetic field strength, is at the heart of many sug- 
gested GRB models (e.g. [6, 22, 9]). Kluzniak and Rud- 
erman, for example, estimate that the magnetic field be- 
comes buoyant at ~ 10 17 G, floats up and breaks through 
the surface of the remnant as a sub-burst. This process, 
winding up the field to buoyancy, subsequent floating up 
and sub-burst, would continue until the energy stored in 
differential rotation is used up (or collapse sets in), ~ 
100 s. Therefore even long bursts, with substructures on 
millisecond time scale set by the motion of the magne- 
tized fluid, could result from a neutron star merger. Once 
the differential rotation has been damped out, the system 
could still continue as a gamma ray burster using the en- 
ergy stored in rigid rotation [21] in case it remains stable, 
or, in case of collapse, one would be left with the 'clas- 
sic' GRB-engine, a black hole plus a debris torus. 
The annihilation of neutrino anti-neutrino pairs into 
electron-positron pairs above the poles of the merger 
remnant has been suggested [7] as a process to produce 
a fireball from thermal energy stored in the disk. We 
find that the main neutrino-emitting region is the inner, 
shock- and shear-heated region of the torus around the 
central object. The prevailing densities lie between 10 10 
und 10 12 gcm~ 3 and temperatures are ~ 2-3 MeV. In 
our simulations we find typical neutrino luminosities of 
~ 10 53 ergs with mean energies of ~ 10, ~ 15 and ~ 20 
MeV (see Fig. 2) for an initially corotating system, and 
slightly higher values for irrotational systems [16]. The 
dominant emission process is the capture of positrons 
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We have further analyzed the neutrino signal expected 
from the event. The neutrino emission is dominated by 
anti-neutrinos produced in positron captures, the total lu- 
minosities lie around ~ 10 53 erg/s and are substantially 
lower than those found in previous investigations. Pre- 
liminary analysis of the neutrino anti-neutrino pair an- 
nihilation efficiency suggests that it is difficult to power 
energetic bursts with energies ^> 10 49 ergs via this mech- 
anism. Whether this energy is enough to power a short 
GRB will have to be clarified by future afterglow obser- 
vations. 



ACKNOWLEDGMENTS 



FIGURE 2. Total neutrino luminosities and mean energies 
(in MeV) for the v e , v e and v x emitted in the coalescence of 
a corotating neutron star binary system of twice 1.4 Mq. The 
abscissa gives the time in ms. 



The calculations reported above were performed on the 
United Kingdom Astrophysical Fluids Facility (UKAFF) 
and on the supercomputer of the Math Modelling Center 
of the University of Leicester. 



onto neutrons which occurs in the hot, neutron-rich in- 
ner regions of the debris torus. The found luminosities 
are lower than those found previously by Ruffert and 
Janka [17] by roughly a factor of 3 to 4. The annihila- 
tion efficiency is currently being investigated in detail, 
but the preliminary results indicate that it will be difficult 
to power Gamma Ray Bursts largely in excess of ~ 10 49 
ergs. 



SUMMARY 

We have performed high resolution calculations of neu- 
tron star coalescences with more than 10 6 SPH-particles, 
using a new, nuclear equation of state and accounting 
for the effects of cooling and compositional changes 
from neutrino reactions by a detailed multi-flavour neu- 
trino treatment. For the most realistic case with negligi- 
ble neutron star spins, we do not see any sign of a col- 
lapse to a black hole and therefore argue that the out- 
come of the coalescence is a (at least temporarily) sta- 
bilized, supermassive, hot, differentially rotating object 
with huge magnetic fields. Such an object could produce 
a GRB in various ways, e.g. via a relativistic electron- 
positron wind (Usov 1994) or via the so-called DROCO- 
mechanism (differentially rotating compact object) sug- 
gested by Kluzniak and Ruderman (1998) where the field 
inside the differentially rotating central object becomes 
locally wound up until it becomes buoyant at ~ 10 17 G, 
floats up and breaks through the surface as a sub-burst. 
This self-limited process would continue until the differ- 
ential rotational energy is used up and would therefore 
continue for several seconds with substructures given by 
the fluid instabilities on a millisecond time scale. 
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